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Abstract 

We analyze the production of triple pseudoscalar Higgs bosons in the decay 
channel of Z — > AAA for light pseudoscalar bosons when the corresponding scalar 
boson is too heavy to be produced by Z decay. Analytic results are obtained both 
at the tree level and at the one-loop level. The branching fraction can be as large 
as 10 -5 which should be detectable at LEP. 



PACS numbers: 14.80c , 13.38 , 12.60.F 



There is essentially no stringent and model-independent limit on the mass of a pseu- 
doscalar Higgs boson, generically denoted by A in this paper. Such a pseudoscalar boson 
always exists in the extended Higgs sector beyond the Standard Model (SM). An identical 
pair of pseudoscalar bosons cannot be produced in pair in the Z decay, as it is forbidden 
by Bose statistics. The potential limit on the mass of a pseudoscalar Higgs boson comes 
from LEP experiments. However, in all the analyses^, 0, 0], the pseudoscalar bosons are 
assumed to be produced by the decay of a physical scalar boson h. For the case when 
the scalar boson is heavy (such as > rriz), no limit on itla has been extracted yet. 
If the scalar partner h is heavy enough, the mode Z — > HA will not be allowed by kine- 
matics. Nevertheless, the channel Z — > AAZ* — > AAl + l~~ is allowed if A is light enough. 
However, its branching fraction was shown to be typically about 10 -8 , too small to be 
detectable for LEP. A pseudoscalar Higgs boson lighter than a b quark can be ruled out 
by b — > Sv40, however the conclusion will be very much model dependent. (Therefore, 
it is still worthwhile to make a direct search at LEP even if the pseudoscalar mass is in 
this light range.) In any case, for a pseudoscalar boson whose mass is heavier than the b 
quark and whose companion scalar boson is too heavy for the decay Z — > HA, the current 
model independent bound on its mass is very weak. 

In this note, we look into another potential discovery channel Z — > AAA for the 
pseudoscalar boson which may be detectable among the rare Z decays. The channel is 
particularly interesting when the lightest scalar is heavier than the lightest pseudoscalar 
boson which can also be an axion. Experimentally, the AAA final states was searched for 
by LEP detectors |TJ assuming that two A's are decay product of a physical scalar boson 
h. A lot of that analysis can probably be borrowed immediately to the case when h is 
of-shell. 

In the popular Minimal Supersymmetric Standard Model (MSSM), at tree level the 
Higgs masses obey relations^ m h < m z < m H , m h < < m H and m H ± > mw, 
where rrih is defined to be the lighter one of the two scalar bosons. These relationships 
are modified when one loop corrections, due to top quark, are taken into account, rrih 
in this case no longer has to be lighter than mz- However, it is still constrained to be 
lighter than about 140 GeV (for tan/5 > 1)0. With the radiative corrections, it is also 
possible |3D that uia < 2 m h- 111 this sense, our analysis is also very much relevant to the 
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MSSM in addition to the more general models. 

The one-loop amplitude for Z — > AAA via the virtual top-quark was roughly esti- 
mated by Li[[7l]. Here we study in details both the tree-level process due to a virtual scalar 
Higgs boson, and the one-loop process due to the top-quark loop. 

In Fig. 1, we illustrate one of the Feynman diagrams that the triple pseudoscalar 
decay mode Z — > AAA occurs through the gauge vertex Z — > Ah*, followed by h* — > 
v4A[[|. Phenomeno logically, one can describe the interaction among the scalar and the 
pseudoscalar Higgs bosons by an effective Lagrangian, 

C = X(V)hAA . (1) 

The coefficient X(V) is related to the vacuum expectation value of the Higgs field, the 
quartic bosonic couplings and also some mixing angles. Its value is about the scale of the 
electroweak interaction. 

The amplitude for Z(p z ,£z) — > A(pi) + A(p 2 ) + A(p 3 ) can be written in term of the 
form factors as, 



M 

where 



We simplified our picture by assuming that the contribution of the lightest scalar Higgs 
boson dominates. The analysis is parameterized model independently such that it would 
be straight forward to adapt our study to a specific model. 

From Eq. (|2|), we obtain the spin-summed amplitude squared as follows, 



x w (l - x w )M 2 z 



2xy -4(1 — z- a) z 2 - 4a 



(1 — x + a — h)(l — y + a — h) (1 — z + a — h) 2 



+ 2 other permutations of x,y,z . (4) 

Here a = M 2 A /M 2 , h = M 2 /M 2 and x = 2 Pl ■ p z /M 2 , y = 2p 2 - p z /M 2 , z = 2p 3 ■ p z /M 2 , 
with x + y + z = 2. The allowed x and y ranges are 

2\fa < x < 1 — 3a , 
l-l( x + d)< y <l-^(x-d), (5) 
withd = (x 2 -4a)5[l - 4a/(l - x + a)]? . (6) 
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The partial width for the channel Z — > AAA is 

In Fig. 2, we demonstrate the branching fractions of the process Z — ► AAA for different 
scenarios, (a) M h = 90 GeV (dashed), (b) M h = 100 GeV (solid), and (c) M h = 110 
GeV (long dashed), for the case X(V) = 100 GeV. We find that the size of the branching 
fraction can be as large as 10~ 5 for lighter h. 

Next, we look at the induced amplitude at the one-loop level. This is potentially 
significant when is very large so that the tree- level amplitude is not important. It 
is also interesting because it only depends on the coupling of the pseudoscalar boson 
with the top quark, and independent of the details of the Higgs self- couplings. We shall 
parameterize the Yukawa coupling of the top quark to the light pseudoscalar boson (A) 
in the following model independent form, 

^Yukawa = Q^l^A . (8) 

The coupling between the top quark and the Z gauge boson is given in the Standard 
Model, 

t e 

® z 4 sin 6w cos 9w 

For the top-induced amplitude of the process Z(pz) — > A(pi)A(p2)A(p 3 ), there are six 
Feynman diagrams. Under charge conjugation, they pair up into three sets, 

P123 = 7Pi -Spa- \{u + l - 13 - 5)p z , (10) 

^i23=Pi23- 1 4 M 2 z ~( 1 + 5)m 2 A + m 2 t +(5p 3 + 1 p 1 )-p z , (11) 

Mvx = VI{M\ - Am 2 - 8/4 3 + Ap 2 123 ) + 4[(p 3 - pxfpz - 2p v lP?J + 2^pJ ■ p 123 . (12) 

The Feynman parameters satisfy a + /3 + 7 + d> = l and < a, [3, j, 5 < 1. In Fig. 3, 
we carefully show the choice of the momentum flow and the corresponding Feynman 
parameters. Our results can be easily produced following such convention. A color factor 
3 has been explicitly included in Eq.(||). It is straightforward to generalize expressions 
(10-12) to other cases 132 and 213 by permutations. The charge conjugated diagrams 
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give equal contributions as one can easily check that A/123 = A/321 etc. To arrange this 
one-loop amplitude in a similar form as (§), we introduce the form factors F t (pi,pj) in 
parallel with F h , 



FV,P3) = 3^||p/ / /Sldadpdj 



(4p z + 8pi) • P132 (4p z + 8p 3 ) • P312 

(/ i 132) 2 (/ i 312) 2 

M|-4mg-8^ 23 + 4pf 23 



(13) 

The above formulas are ready for numerical integrations. However, for the purpose of 
illustration we only extract the leading contribution in the large mt limit even though the 
correction can be of order of M^/mf « 0.25. We also remove irrelevant constant terms 
from F f . Such constant terms cannot contribute to the overall amplitude for a physical 
polarization e of the Z boson. We obtain, 

■ (14) 

Unfortunately, such a top-loop induced amplitude is so small that it produces, by itself, a 
negligible branching fraction for Z — > AAA below 10~ 10 , even we assume a SM coupling 
gt = {\/2Gp)^rat. This is much smaller than the previous rough estimate by many 
orders of magnitude. More likely, the signal of Z — > AAA comes from the Higgs mediated 
process. 

Since one requires the scalar boson h to be light enough (such as 90 GeV) in order to 
get a large branching ratio, one may also consider the alternative production of e + e~ — > 
Z* — > HA —* AAA. This possibility is already covered in some of the Higgs search 
analysis 0, §. 

As the accumulated events of Z — > hadrons among the four LEP groups have reached 
10 7 , A branching ratio of 10~ 5 is potentially detectable. The main difficulty seems to be 
finding a clear signal with high efficiency for such events. If the pseudoscalar boson is 
heavier than 2m& ~ 10 GeV, then presumably it will decay dominantly into six b quarks. 
In MSSM[PJ, for bosons decay into bb about 90% of the time and about 6 — 8% into t + t~. 

For the case 10 GeV ~ 2m^ > > 2m T w 3.5 GeV, the A boson can decay 
dominantly into six r leptons or six charm quarks. The two modes are competitive 
with each other. One can search for t + t~ plus four jets or t + t~t + t~ plus two jets 
or t + t~t + t~t + t~ in increasing detection efficiencies. The answer will depend on the 
relative fraction between t + t~ and cc final states. 
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For the channel Z — ► AAA — ► bbbbbb, the clear signal can be a number of 6-tagged jets. 
Similar signal was searched before in the previous Higgs search analysis [|IJ for lighter on- 
shell scalar boson h using the same AAA final state. It was concluded]!] that the current 
limit of this branching ratio is at about 10 -4 level. With more recent data, this limit may 
be improved by a factor of 3 or more with improved statistics. To improve this further, 
one probably has to increase the efficiency in the identification of six jets from the three 
A bosons and the efficiency in b-tagging. Typically a prize (of about 20 — 30%) has to 
be paid to impose tight cut to reject 3, 4 or 5 jets events. The efficiency will be higher 
for lighter pseudoscalar. In addition, one has to pay a prize for b-tagging. The current 
b-tagging efficiency of LEP detectors is roughly about 20% per jet. Even if one tags only 
3 out of 6 jets, the prize is already quite severe (20 x (0.2) 3 = 16%). These two effects 
combine to give 3 — 5% efficiency of identifying Z — > AAA. (OPAL|l]] quoted 6 — 11%, 
but with rather high background). 

In a general multi-doublet extension of Standard Model, it is possible that the pseu- 
doscalar decays into r leptons or b quarks with similar branching fraction, in that case, 
the best modes to discover pseudoscalar boson may be bbr + T~T + T~ or bbbbr + r~ final 
states. H As far as we know, these modes have not been seriously searched for by LEP 
yet. 

Our analysis indicates that there is a good chance that one can detect signal of pseu- 
doscalar boson in Z decay with a branching ratio of about 10~ 5 . It is certainly far from 
trivial; however, the reward is that the pseudoscalar boson may be already there in the 
data waiting to be uncovered. 
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Figure Captions 

1. One of the tree-level Feynman diagrams for the process Z — > AAA via the scalar h 
Higgs boson. Two other diagrams are obtained by permuting the momenta. 

2. Predicted branching fractions of the tree-level process Z — > AAA for different sce- 
narios, (a) M h = 90 GeV (dashed), (b) M h = 100 GeV (solid), and (c) M h — 110 
GeV (dot-dashed). We have set (V) = 100 GeV. 

3. One of the one-loop Feynman diagrams for the process Z — > AAA via the virtual 
top-quark. The choice of momentum flow and the Feynman parameters are clearly 
labeled. Two other diagrams are obtained by permuting the momenta. 



A 



\ 



\ 



Fig. 1 One of the tree-level Feynman diagrams for the process Z — > AAA via the scalar h 
Higgs boson. Two other diagrams are obtained by permuting the momenta. 
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3 One of the one-loop Feynman diagrams for the process Z — > AAA via the virtual 
top-quark. The choice of momentum flow and the Feynman parameters are clearly 
labeled. Other diagrams are obtained by permuting the momenta. 



